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Abstract. - We show that the distribution P(g) of conductances g of a quasi one dimensional 
wire has non-analytic behavior in the insulating region, leading to a discontinuous derivative in 
the distribution near g = 1. We give analytic expressions for the full distribution and extract 
an approximate scaling behavior valid for different strengths of disorder close to g = 1. 



It is well known that the probability distribution P(g) of the (dimensionless) conductances 
g of a disordered conductor is Gaussian in the deeply metallic regime and log-normal in the 
deeply insulating regime Q . It has been proposed recently Q that for a quasi one dimensional 
(ID) wire with mean dimensionless conductance g <C 1, P(g) has a sharp cut-off beyond .9 = 1. 
To be sure, in one dimension, when g < 1, P(g) drops discontinuously to zero at g = 1. In 
higher dimensions the discontinuity will be smeared in some way, but it is difficult to say how 
much of it will survive. Using the saddle point method developed in ||, we are able to show 
here that nonanalytic behavior remains near g = 1 at least in quasi-lD systems. One should 
note that the exact results for the mean and variance of the conductance, obtained within the 
non-linear sigma model [jj), do not give any clue on the abrupt change of P{g) near g = 1. 
The existence of the sharp cutoff for all g -C 1 has recently been confirmed numerically, while 
the exact shape of the distribution very close to g = 1 seems to be different from a log-normal 
distribution 

In the present work, we focus on the insulating region in quasi ID near \g — 1| <C 1, in 
order to understand better the nature of the unexpected sharp feature in the distribution. For 
simplicity, we will restrict our discussions to the unitary case of broken time reversal symmetry. 
We propose from a simple generalization of ^ that for large enough disorder the distribution 
has a non analyticity near g — 1, giving rise to large discontinuities in its derivatives. We 
obtain e.g. that in the insulating limit, the leading contribution to the distribution P(g) at 
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1.9 - 1| < 1, is given by 



(yf{v - i)) - $ (yf{v x - i)) + C; .g > 1 



p (flO «S>\ / \ (!) 

' *(^))-#(Vr(i*-£))+C7; 5<l + «, 

where $ is the error function, v = cosh -1 y/l/(g— 1) , i/i = cosh -1 y/2/g, and a= 1/ cosh 2 (3/2r) 

The disorder parameter T = (,/L, where £ is the localization length and L is the system 
size (r> 1 would correspond to metals and r< 1 corresponds to insulators). The term C is 
independent of g, but depends on T. The discontinuity is at g = 1 + a. For the case considered 
here, r < 1 and a ~ e~ 3 / r . From (Q) it follows that P' = dP(g)/dg has a discontinuity at 
3 = 1 + 0!, with P' ~ — e 2 / r very large for g > 1 + a (growing exponentially with increasing 
disorder) and \P'\ <C 1 for g < 1 + a. Note that in the limit r — > 0, the distribution would 
have an essential singularity at g = 1. 

The above results are obtained assuming r « 1, describing the insulating limit. As T 
increases, but still in the insulating regime, the dominant contribution remains the same with 
renormalized parameters, and our results remain qualitatively valid. In the metallic limit the 
singularity is absent, giving rise to the possibility that the singularity disappears across the 
crossover region (r ~ 1/2), where our current approximations are not valid. 

In addition to the sharp structure, our results suggest an approximate scaling behavior for 
different strengths of disorder near <? ^ 1. As is evident from ([!]), in the expression for the 
ratio P(g)/P(l), the dominant dependence on the disorder parameter T and on g — 1 appears 
only in the combination s/T(y — l/2r). Our numerical results agree with this approximate 
scaling. We also find that the distribution P(g) for g < 1 is a slowly varying function of g, 
which is approximately constant close to g — 1. This leads to an exponential distribution for 
P(ln g) in this region, rather than the log-normal distribution proposed in [Q, and agrees with 
various existing numerical results ||. For Y <C 1, the numerical results suggest that the 
exponential form crosses over to the log-normal form for g 1. 

We now briefly discuss the details of the method which is based on a simple generalization 
of The probability distribution p(X) of the N variables A-i, where Xi are related to the 
transmission eigenvalues Ti of an N-channel quasi Id wire by A, = (1 — Tj)/Tj, satisfy the well 
known DMPK equation p0||, whose solutions in the metallic and insulating regimes can be 



written in the general form 1 1 



p(A) = |exp[-/3fr(A)], (2) 

where Z = J FJ^ dAj exp[— (3H\ is a normalizing factor independent of A^, H{\) may be inter- 
preted as the Hamiltonian function of N classical charges at positions Aj, and (3 = 2 for the 
unitary case. Since the dimensionless conductance is given by 



the probability distribution P(g) can be written as 



i roc j rco N N 



p(\). (4) 
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Following {§, we define a "Free energy" 

^) = £t^+/^ (5) 



N 

IT 



such that the distribution can be written as 

N 



1 f°° c\t f°° 

P(g) = % / Y 6 " 9 / II dA * cx p l~ F W • ( 6 ) 

It has been shown in Q that separating out the lowest eigenvalue Ai and treating the rest 
of the eigenvalues as a continuum beginning at A2 , the distribution can be obtained within a 
generalized saddle point approximation, given by 

/•OO poo 

P{g) = / dAr / dA 2 e- s , (7) 
Jo J\i_ 

where 

' (g-F') 2 + F° (8) 



2F" 

is the saddle point action, obtained from a saddle point free energy 

F sp = F° + (ir)F' + ^F". (9) 

The insulating region in this model is given by Xi 3> 1, and also X2 3> x\, where sinh 2 Xi = \. 
The saddle point free energy terms in this limit were calculated to be 

F°(x u x 2 ) » 2T 2 xl -6Tx 2 +Txj - - ln(xi sinh(2a;i)); (10) 

(11) 



F" = — 



1 



sinh (2x 2 ) 



cosh 2 x\ 
11 1 



3 4x\ sinh 2 (2a;2) 



(12) 



Since 1/2|_F"| is exponentially large in the insulating regime, the saddle point solution is given 
by putting the coefficient of the term in the action S equal to zero, namely g = F' , which 
then gives 

coshicisp = — . (13) 

As pointed out in the saddlepoint solution is valid only for g < 1 since coshxi > 1, while 
the g > 1 region is determined by the boundary values of x\ =0 and x 2 = 2/irT. While this is 
a good approximation on both sides of g = 1, the fact that the boundary of the saddle point 
solution is at g = 1 makes it possible that the approximation is not accurate enough very 
close to g — 1. We will show below that at the next level of approximation suggested in 
the region close to g = 1 can be better described by saddle point solutions valid on both sides. 
This improved approximation immediately leads to the non-analyticity mentioned before. 

It was shown in Jl3| that while separating out the lowest eigenvalue gave qualitatively 
good results both in the metallic and insulating regimes, it is important to separate out one 
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additional eigenvalue to obtain good agreement with numerical and available exact results in 
the insulating and crossover regimes because the separation between the eigenvalues becomes 
large. The extension is straightforward, and the distribution P(g) now has an additional 
integral over A3, which is the beginning of the continuum that represents the rest of the 
eigenvalues. It is clear that the insulating limit is characterized by X3 ^> X2 S> 1 for typical 
values of x% , and that F' in this limit will now be given by [ p"3| 

F' « \— + \—. (14) 

cosh X\ cosh X2 

F" now has the same expression as (12), with X2 replaced by X3, while F° is given by ( [To| ) 
with X\,X2 replaced by X2, X3 plus additional terms involving x\ 

F°(x 1 ,x 2 ,x 3 ) = F°(x 2 ,x 3 ) + Yx\ - - ln(xi sinh(2ari)). (15) 

The saddle point solutions for 0:3 and X2 are now given by X3 Sp = 3/2T, which is indepen- 
dent of g and 

coshx 2sp = , (16) 

J g — 1/ cosh 2 x\ 

obtained again from g — F' . Since X2 3> 1, the additional term in F' is exponentially small, 
and usually negligible. However, close to g — 1, this is the term that allows the saddle point 
solution to be valid on both sides of g = 1. The fluctuation correction to the X2 integral 
contributes an additional factor e 2x2sp , and we are left with the final integral 

Pig) oc f ^ dx 1 e- rx ' +lnxi ' rx ^p +3x2s ". (17) 

We have used finite limits in the integral because the saddle point solution of X2 must be 
real and also we require that x 3sp > X2 Sp - This gives x\ m i n = i?e[cosh~ 1 \J\j{g— a)], a 
I cosh 2 (3/2r), and xi max = cosh -1 \/2/g. The crucial point is that the restriction x\ m in is 
zero for all g > 1 + a, but becomes finite for g < 1 + a, with a discontinuous slope. It is the 
universal nature of this constraint which will eventually lead to the sharp changes in P(g) at 
g=l + a. 

In the limit g — ► 1 and r 1, the upper limit of the integral x\ max <C 1. We can therefore 
neglect the x\ term in (JT^). Changing variables to x 2 and rewriting lnxi, we obtain a simple 
integral 

P{g) oc [ 2maX dx 2 e- rx * +X2 (18) 

J V\ 

with v\ as defined after ([!]), and X2max = v for g > 1 + ot, but 3/2r for g < 1 + a. The 
integral immediately leads to Eq. (|l|). It is easy to see from ([j]) that the first derivative has 
a discontinuity at g = 1 + a and is large, of order e 2 / r for g > 1 + a. Figure [j] shows the 
sharp changes at g = 1 + a according to Eq. (|l|). Since v\ in this region is only weakly 
dependent on g, this also leads to an approximate scaling behavior of P(g) as a function of 
z = y/T(v - 1/2T). 

While the g-dependent contribution comes from the X2 integral where X2 S p < %3spi there 
is a contribution to P(g) from integral over X2 > xs sp . The X3 integral in this case is given 
by the boundary value at X3 = X2, and the fluctuation correction to the X2 integral becomes 
of order unity. The resulting integral leads to a g independent but T dependent term giving 
C«e 3 / r [l-$(2/v / T)]. 
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Fig. 1 - P(g)/P(l) from Eq. (Q) for different values of I\ 



In the above, we considered the case r<l, where x 3 x 2 S> 1 and only the leading order 
terms were kept in the free energy. For weaker disorder but still in the insulating side, the next 
leading order terms will become important. First of all, the In £2 term arising from F°(x2, X3) 
in ( |l5| ) will start to contribute. Second, the next leading order contributions to F° will 
contribute some more x 2 and In £2 terms. In general, as long as the limit 13 > 12 > 1 remains 
valid, and the pcrturbative expansion of the free energy remains a good approximation, the 
integral ( |l8| ) can be generalized to be of the form 

P{g)oc [ X2m " X dx 2 x p 2 e- arx i +bX2 (19) 

J V\ 

where a, b and p are parameters which in principle can be calculated within our framework 
by keeping higher order terms in the free energy. For p — we get back ([!]) with renormalized 
parameters. For finite integer p the integral can again be evaluated simply in terms of the 
error function. For non-integer values of p, the integral is more complicated but it is clear 
that the qualitative features of the results will remain the same. We expect that our results 
will remain qualitatively valid near the crossover region on the insulating side (T < 1/2). 

In order to check some of our approximate expressions, we have calculated P(g) numerically 
for a quasi ID wire. Recent techniques allow us to obtain reasonable statistics for N = 5. 
Unfortunately, it was not possible to obtain enough data for either r < 1 or 1 < g < 1.01. 
Nevertheless, ( |l9|) should be approximately valid in the region numerically accessible to us 
and the data can be used to check if the approximate scaling expected from (|l9| ) holds 
qualitatively. In Fig. 0, circles and squares show the numerical results for different values of 
average conductance g. We use p — 2, T e = aT and b as fitting parameters to fit the data 
for different g using (19) (for our range of parameters, the constant C of (0) turns out to be 



negligible). Agreement with the data is good for g — 1 < 0.15. In Fig. || the same data are 
plotted as a function of the scaling variable s = exp[— \/Tl(v — 6/2r e ] shifted by s(l + a), 
taking into account the fact that the singularities occur at g = 1 + a(T e ). Again the scaling 
is quite good. We have checked that the scaling fails for larger values of r e , which is expected 
because it is close to the crossover regime (r ~ 1/2) where our perturbative expansion of the 
free energy starts to become invalid. 

A non-analyticity in P{g) was obtained in 0] within Random Matrix Theory (RMT), based 
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Fig. 2 - Numerical data for N = 5 with different disorder compared with theoretical predictions 
using (p"s|). Data for g ~ 0.3 were fit (solid line) with T e — aF — 0.292 and b — 2; g w 0.4 were fit 
(dot-dashed line) with T e = 0.346 and 6 = 1.8; g ~ 0.5 were fit (long dashed line) with r e = 0.391 
and b = 1.6. In all plots, p — 2. 



on the assumption that contributions from only the smallest two eigenvalues are important 
in the insulating regime. All derivatives of P(g) in such a model diverge at g = 1. However, 
a two-eigenvalue calculation within DMPK (obtained by neglecting the continuum in our 
model) gives rise to an essential singularity at g = 1, so that there is no discontinuity in 
the derivatives at that point. Therefore, it is not simple to see if the trivial discontinuity 
in the derivative mentioned in the introduction will persist beyond ID. We emphasize that 
our present results for quasi ID are fundamentally different from both of the above models: 
whereas all derivatives of P(g) diverge for the RMT model and vanish for the two-eigenvalue 
DMPK model, they are finite for our quasi ID system. Also, in contrast to the above cases, 
the singularity for the quasi ID system is not exactly at g = 1, but is shifted as a function of 
disorder. Finally, the RMT results do not satisfy the approximate scaling of Fig. |[ 

We note that there is no phase transition in quasi ID. It is therefore important to ask if 
a non-analyticity in P(g) in the absence of a true phase transition violates any fundamental 
principle. Our framework can be considered as an electrostatic problem in one dimension 
involving charges with repulsive interactions and confinement potentials. Our free energy 
and the density of charges are all analytic, ensuring that e.g. the total energy will be an 
analytic function of disorder. However, the conductance is a complicated function of the charge 
distribution, and the non-analyticity appears only in the distribution of the conductances 
where the charges have to satisfy certain constraints. There is no restriction on the analyticity 
of such a distribution. 

In summary, we have shown that the distribution of conductances in quasi ID systems 
in the insulating regime has a non-analytic behavior near g = 1. The non-analyticity gives 
rise to very sharp structures close to g — 1, with finite discontinuities in its derivatives. The 
presence of similar structures in higher dimensions in numerically studied systems |?], || 
gives rise to the possibility that such non-analyticity might be present in the conductance 
distribution in higher dimensions as well, having important consequences for the Anderson 
transition. Whether the non-analyticity discussed here disappears abruptly at some critical 
value of r, or smoothly as T is increased beyond 1/2 remains to be clarified. It is conceivable 
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Fig. 3 - Appropriately scaled data from Fig. ^ showing approximate scaling as a function of s = 
exp[— y/T2(v — b/2T e ], shifted by s(l + a). The theoretical lines from Eq. (19) for different T e (solid, 
dot-dashed and long-dashed lines for T e = 0.292, 0.346 and 0.391, respectively) do not fall exactly on 
each other, showing the approximate nature of the scaling. 



that an abrupt disappearance of the non-analyticity is a signature of the Anderson transition. 
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